Glass phases of flux lattices in layered superconductors 
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' We study a flux lattice which is parallel to superconducting layers, allowing for dislocations 

, and for disorder of both short wavelength and long wavelength. We find that the long wavelength 

' disorder of strength A has a significant effect on the phase diagram - it produces a first order 

t-H \ transition within the Bragg glass phase and leads to melting at strong A. This then allows a Friedel 

» ! . scenario of 2D superconductivity. 
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■ The influence of disorder an the behavior of flux arrays in superconductors is intensively studied in recent years. 
t-H [ A model suggested by LarkinEJ in which random, forces act independently on each vortex was applied to conventional 

superconductors. Usins,.tbe arguments of Refill itpjwas shown that the flux array is collectively pinneda, forming a 
' ' ' vortex glass ( VG) phaseatffi. However, experimentsQB on weakly disordered samples, reiveaUong-range order of the flux 
array, much beyond the characteristic Larkin length. This phenomena is accounted fom'E30 by an elastic dislocation- 
free theory of the flux lattice in weak random potential which predicts an algebraic decay .nf. the translational order 
and the existence of divergent Bragg peaks in the glass phase (Bragg glass). It was arguedlJO that the Bragg glass 
<*S !■ phase is stable against formation of dislocations in a finite range of the phase diagram. __ _ ___ 

The phenomena of melting of the flux lattice is of considerable interest. Melting has been observedLn3'o : E^loll3 as 
either a transition into a flux liquid phase or into a glass phase with a higher critical current. A model for melting, 
allowing for both disorder and dislocationsErEj was recently studied. The model considers flux lines parallel jta 

Sand confined between superconducting layers and allows for dislocations. This model was studied without disorderlill 
leading to flux melting at a critical tempera ture . T c which is about factor ~ 2 from the solution of the more fundamental 
r T^ \ system in terms of superconducting phasesHLa; the latter allows also for flux loops and overhangs. Since disorder has 
drastic effects on melting we expect that the model in terms ofjkrx displacement is a reasonable starting point. 
The solution of Carpentier, Le-Doussal and Giamarchi (CLG)E3 has shown explicitly that short wavelength disorder 
i^i - combined with dislocations leads to melting at a finite value of the disorder strength. CLG havej-used Replica 
Symmetry Breaking (RSB) methods as well as Renormalization Group (RG). They have also showntSEj that this 
04 melting is compatible with a Lindemann criterion. 

In the present work we allow for an additional term in the CGL model. This term is generated by RG within the 
CLG model and leads to significant effects on the phase diagram. Using RSB methods, we show that long wavelength 
disorder of strength A leads to a first order transition within the Bragg glass phase and as A increases it leads to 
| melting. We find that the A induced melting is inconsistent with a universal Lindemann criterion. Finally we consider 

■ the quest for the Friedel scenarioEj in which a layered superconductor becomes a set of decoupled two-dimensional 
' (2D) .superconductors. This scenario fails in pure supercondutorsErtJ, but is possible with some constraints in parallel 

s * i fieldso and in special modelscj. With disorder which affects interlayer coupling the Friedel scenario becomes feasible 
in presence of a rnpjted flux array. 

The modeOEjO consists of layers with interlayer spacing I where modulation in the flux line density couples to a 
random potential. We consider a Hamiltonian with two types of random potentials, 
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O ■ -2Re(Q(r)e^)} (1) 



H = j d 2 r^[^(V$ l (r)) 2 -7 ?l (r)V$ l (r)- M cos($ J (r)-$ 4+1 (r)) 



with Gaussian disorder correlations < Q(r)Q(r') >= 4Tg5i.jS(r — r') and < r]i(r)rjj(r') >= TASij5(r — r') where T 
is the temperature. Here $j(r) stands for in-plane displacement of the vortex line in the i-th layer, c is an in-plane 
elastic constant, g measures disorder with Fourier component w 2it /a where a is the flux periodicity parallel to the 
layers, while A measures long wavelength disordertj. The /i term is the coupling between layers which allows for 
dislocations. In the pure systemlill thermal fluctuations lead to melting, i.e. /i is renormalized to zero, at T c — Aire. 

Dimensional Imry-Ma arguments are useful to check the stability of an ordered phase which is a d dimensional 
elastic medium. In a domain of size L the elastic energy is oc L d ~ 2 , the short wavelength disorder (after averaging 
the square) is oc L d / 2 while the long wavelength disorder is oc £( d_2 )/ 2 . Thus short wavelength disorder is relevant at 
d < 4, while the long wavelength disorder is marginal only at d = 2, i.e. the latter is consistent with long range order 
in d=3. 

To average over disorder we start with the replicated version of Hamiltonian (Eq.(l)) which includes all relevant 
terms generated by renormalization, 
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H = J rfM|E[( v$ "( r )) 2 -^ cos W( r )- $ i+i( r ))] 

i.a 

^[|v$»(r)V$f(r) - 1C os(^(r) ~ *? + i(r) - *J(r) + $ 4 b +1 (r)) 

- fl cos($?(r)-$t(r))]} (2) 

where a = l...rt is the replica index. Note in particular the 7 term which was not considered by GLC; this term is 
generated in second order RG from the p term. Since it couples different replicas it can lead to RSB, i.e. this term 
leads to distinct phenomena and should be included in the full Hamiltonian. 
We consider the variational free energy F var = Fq+ < H — Ho > with 



(3) 



where the Greens' function G~ b (q, q z ) is determined by an extremum condition of F var and q, q z are Fourier variables 
for r and i, respectively. 

Defining the inverse Green's function in the form G~ b (q, q z ) = 5 a bG ~ 1 (q,q z ) — <j a b — Aq 2 with ^2 a a o-b — and 
a a b = 2(1 — cos q z )a2 b + <J 9 ab we obtain the the self-consistent equations in the form 

G Q 1 (q,q z ) = cq 2 + 2p(l-cosq z ) (4) 

* 9 ab = 2 9 eM-~B 9 J (5) 
a2 b = 7 exp[-B2 b ] (6) 

p = fiexp[--B aa ] (7) 

where we define 

B 9 ab = 2Tj2lGaa(q,q z )-G ab (q,q z )} (8) 

Kb = 27^(1 -cosq z )[G aa (q,q z ) - G ab (q,q z )} (9) 

B aa = 2T^(1 - cos q z )Gaa(q,q z ) ■ (10) 

Here p is the renormalized coupling between layers which is determined by the diagonal B aa ; p = signals a 2D 
phase, i.e. correlations in the z direction are lost and the flux lattice has melted. .—. 

We study the full RSB solution of the saddle point equations (5-7). The method of RSBcJ employs a representation 
of hierarchical matrices such as <j' b in term of functions a gn (u), and similarly B 9 ' b is represented by B gn {u) with 

< U < 1. r-| 

We define two order parameters for RSB, [c gi7 ](it) = ucr gn (u) — J Q l a gn (v)dv. Using the inversion formulae^ for 
G a b and integrating over momenta we obtain 



]-B g<1 (u) = \g gtl {u) - [ %g g ,~f(s) (11) 

Z 11 J „. s 



g g {u) = -T\n[-{p + l + w)] 

g^^ggi^-fip+l + w)- 1 (12) 

Here p{u) = [w(u)(w(u) + 2)] 1 / 2 , w(u) = ([a g ](u))/2m(u), m{u) = p+ [cr 7 ] (u), A is a cutoff of cq 2 (A 3> m(u), [a g ](u)) 
and T = T/T c with T c = Aire. 

By differentiating Eqs.(5-6) we obtain two coupled differential equations for the RSB functions m(u) and w(u) 

2T dm d . mp(l + w + p) . 



u du du p{\ + w + p) + Q(w + p) 

^(Q + w) = ^[ mp ^ + w) } (13) 
u du du Q + p 
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where Q(u) = m(dw / du) / (dm / du) . 

A general solutipa of these equations is rather difficult, so at first we consider special limits. When 7 = we recover 
the CLG solutionis. Within the 3D Bragg glass phase < [$;(r) - $;(0)] 2 >~ In r so that positional correlations 
decay algebraically and long range order is weakly destroyed. The Bragg glass phase undergoes a continuous melting 
transition (for A = 0) at g/ [i = 2/eT as shown in the A = plane of Fig. 1; for A ^ the transition becomes first 
order. Thus, the Bragg glass phase, due to both disorder and dislocations, melts into a 2D phase with p = 0. 

Consider next thecase g = 0, hence w(u) — 0; the solution in this case is formally similar to that of a 2D disordered 
Josephson junctior£3. Eq. (13) yields then (1 — 2T /u)m'[u) = 0, i.e. m(u) is a one step function, with the step at 
u = 2T. Since u < 1 the onset of this solution is at T — 1/2. i.e. at T = T c /2. Eqs. (6,7) determine the jump in 
[<7 7 ](w) from zero (u < 2T) to a value ct 7 at 2T < u < 1, where 



p + (7° ( - 7 ^l/(l-2f) 

— e 



A 



1 (4eTl) (14) 



A+l/2 I A„ml \ ' A 4 



1/2 



A 



(15) 



where A = 4irT 2 A. This solution is valid for A < 1/2 and T < 1/2. For 7 of order \x, near the T = 1/2 transition /2 
is finite while /x + c 7 vanishes. Thus cr" < is finite up to T = 1/2 and vanishes at T > 1/2, i.e. the transition is of 
first order. When A > T within this phase a° changes sign and becomes positive. 

The phase at A < 1/2 and T < 1/2 is a coexistence phase-it has both long range order (p ^ 0) and glass order 
(cr° 7^ 0). (As noted above this is consistent with the Imry-Ma argument). At A = 1/2 we find a disorder driven 
transition where p vanishes continuously, leading to a 2D glass phase at A > 1/2. 

We note also that a replica symmetric solution is possible for 1/2 < T < 1 with 

i =e -( e A + v^y<^ (16) 

i.e p ^ for T+ A < 1, T > 1/2 as shown in Fig. 1. Comparison with Eqs. (14-15) shows that p is also discontinuous 
at the T = 1/2 transition. 

Finally we consider the case where both g and A are finite. We can demonstrate the existence of a first order 
transition at small g by showing a coexistence of two solutions. The first solution is an expansion near the 5 = 
solution, i.e. w(u) <C 1 with [a g ] — O(g), [07] = O(l); for u < T/2 the solution for [a g ] is similar to the 7 = case, 
i.e. [(T g ](u) ~ u 2 for small u, consistent with Bragg glass correlations. This solution is valid (assuming fi ~ 7) if 

1-3A + 2TA 

g/A < (7/A) (1-2T)(1-2A) (17) 

i.e. for weak coupling g/A, 7/A <C 1 this expansion breaks down close to the transitions at T = 1/2 and A = 1/2. The 
second solution is an expansion around the 7 = solution with [er 7 ](u) -C p. This leads to [a g ] — 0(g 2 ), [<7~] = 0(g) 
and is valid for A < T. Thus for small g there is a two solution regime which implies a first order transition at some 
T < 1/2. We indicate this transition by a spaced dashed line in Fig. 1, though we do not know its precise location. 

As shown in Fig. 1, we find that the main feature of the CLG scenario is valid- for small disorder the Bragg glass is 
stable, while at large disorder, which can have either short or long wavelength, dislocations are enhanced by disorder 
and lead to melting. _ 

These analytic results for melting allow us to test the jLipdemann criterion, which is of common useu. For the 7 = 
case, CLG consider a Lindemann criterion of the formt§E3 < [$ i+1 (r) — $i(r)] 2 ji-j= c 2 L , with average done in the 
elastic limit, i.e. the cosine of the fi term in Eq. (1) is expanded, This criterion leadalj to a reasonable value of cl ^ 1. 
For the g = case an elastic limit leads to an expansion of both the /z and 7 terms in Eq. (2) so that RSB is not 
induced. Since long range order is present the Lindemann criterion is < $| (r) >= c 2 L ; however the replica symmetric 
solution yields < $|(r) >= (A + T) ln(A//i), i.e. at melting c 2 L as ln(A//i). since A//x depends on the anisotropy of 
the system the Lindemann number cl is non-universal. __ 

In order to relate the phase diagram to the actual magnetic field B we need to identify c by the elastic constantsoO 
which are dispersive, c| 4 « cf x = (B 2 /4tt)/(1 + \ 2 q 2 + \ 2 q 2 ); here A, A c are penetration lengths and e = A/A c < 1 is 

the anisotropy (cjj 4 has a smaller second term which is neglected here). The behavior near melting is dominated by 
q — > and q z w 1/2. The lattice periodicities satisfy / = ae if I > d for weak fields, i.e. B = <fr /a 2 e < <fi e/d 2 (d is the 
spacing of the superconducting layers, which is the lower bound on the interlayer spacing / of the flux lattice, and (f>o 
is the flux quantum), or I = d for strong fields, B = 4>o/ad > foe/d 2 . By rescaling the x, z coordinates we identify 
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Aire « a^/(47r 2 AA c ) - B~ 1/2 
4ttc w «i^/(47r 2 A 2 ) 



5 < o e/rf 2 



(18) 



T c (= 47rc) is smallest for large fields, but even then its value is < lOOT^T for typical CuOi superconductors only very 
near the superconducting transition where A diverges. Thus it is the disorder induced melting which is relevant to 
experimental data. _ . 

Since A couples to V$i(r)E£l it is B independent so that A = A/4ttc 2 ~ B for small fields and A ^constant for 
strong fields. Thus the A induced melting can be induced by increasing the magnetic field if A = 1/2 is achieved 
for weak fields. On the other hand, for the g induced melting g ~ l/a 2 E9, and by using the Cqq elastic constantH we 
identify /i = a 2 (poBe 3 /{4:Tr 2 l(8ir\) 2 }. For weak fields the melting temperature is T^ 1 = eg/(8TTCfi) ~ B z / 2 while for 
strong fields T," 1 ~ B 3 . In this case melting is induced by increasing B for both weak or strong fields. 

We note finally that any disorder which melts a flux lattice parallel to superconducting layers is a route to the Friedel 
scenariccj. In this scenario the layers remain superconducting, so that the system behaves as a 2D superconductor. 
This scenario is not valid in pure layered aiiperconductoruO, it is. valid for parallel fields with restricted system 
parameters (e.g. large vortex core energy)E3 or in special modelsEi In the present study, disorder which affects 
interlayer coupling leads to melting which decouples the layers; assuming weak intralayer disorder the layers remain 
superconducting. Thus we have a model system which exhibits the Friedel scenario. 

In conclusion we have shown that a new interaction term, generated by RG, leads to a significant role of the 
long wavelength disorder. This interaction extends the CLG results to the more complex phase diagram of Fig.l. 
We find that the Bragg glass is stable for weak disorder of either short or long wavelength. The long wavelength 
disorder induces a first order transition within the Bragg glass phase; it also leads to melting which is inconsistent 
with a Lindemann criterion. We propose that experiments with parallel magnetic fields can test the present theory 
of melting as well as test the possibility of 2D superconductivity. 



Acknowledgments: We are grateful to E. Zeldov, Y. Y. Goldschmit, A. Kapitulnik, M. V. Feigel'man and D. R. 
Nelson for useful discussions and to T. Giamarchi, T. Nattermann, P. Le Doussal and V. M. Vinokur for discussions of 
their works and for illuminating comments. This research was supported by a grant from the Israel Science Foundation. 



4 



Figure Caption 

Phase diagram of layered flux lattices with disorder of short wavelength (with strength g) and of long wavelength 
(with strength A). The dashed line is the line of first order transition (discontinuity of the interlayer coupling and of 
the glass order); the spaced dashed line is its approximate extension to g / 0. The 3D glass regime at g = has long 
range order while the 3D Bragg glass at g =^ has algebraically decaying positional order. 
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